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Abstract 


A diffusion-limited transient mathematical model for a lithium/air cell, with the air cathode pores flooded with an organic electrolyte, has been 
developed. During cell discharge, the cathode pore radius profile is reflective of the distribution of the lithium peroxide product in the cathode. The 
cathode pore radius profile has been predicted as a function of time, current density, oxygen gas pressure, and cathode thickness for an assumed 
initial porosity and average cathode pore size. Transient concentration profiles of the dissolved oxygen in the electrolyte were also determined. 
Capacities of the lithium/air cell were predicted and compared favorably with literature experimental results. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Background 


Metal/air batteries are unique in that the air cathode active 
material, oxygen, is not stored internal to the cell/battery system. 
Oxygen is provided by access of air to the porous carbon cathode 
electrode for the half-cell oxygen reduction process. Table 1 
shows the characteristics of some metal/oxygen battery couples 
as given by Abraham and Jiang [1]. 

From Table | it is obvious that the lithium/oxygen battery has 
the highest theoretical specific energy. However, as the battery 
discharge process progresses, the metal/air batteries get heav- 
ier with time, thus the true theoretical specific energy for this 
type of battery for the entire discharge period is the average of 
the values of the “excluding O2” and that of “including O2”. 
For lithium/oxygen, this average theoretical energy density is 
8170 Whkg~!. This value is higher than the heat of combustion 
lower heating value (LHV) energy density of 5524 Wh kg”! 
for methanol/air, but, less than the LHV energy density of 
11,860 Wh kg"! for gasoline/air [2]. Thus, with metal/air batter- 
ies typically achieving 33% of their theoretical specific energy 
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values or higher, this type of battery has the potential to be the 
lightest power source for powering portable electronic equip- 
ment, unmanned aerial vehicles, camping equipment, or any 
equipment where air is present [3]. 

The lithium/oxygen organic electrolyte battery differs from 
the aqueous metal/air batteries in that the oxygen reduction prod- 
uct is insoluble in the organic polymer electrolytes [1,4,5]. Work 
performed using a rechargeable lithium/oxygen battery, using a 
lithium ion conducting organic polymer electrolyte, led to the 
conclusion that the primary cell reaction is given by the follow- 
ing equation [1]: 


2Lis) + Org) = Li2Or,s) (1) 


From the oxygen consumption measurements and experimental 
data analysis [4], the reaction given by Eq. (2) was found to also 
occur during the lithium/air cell discharge for certain organic 
electrolytes and at higher discharge rates 


2Lig) + 40x) = Li20%) (2) 


Read et al. [4,5] have shown that for a given flooded-electrolyte 
design, the organic electrolyte formulation and the oxygen con- 
centration have the largest effects on the battery discharge 
capacity and rate capability. It was concluded that the cells 
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Table 1 
Characteristics of some metal/oxygen battery couples 


Metal/O» couple Ideal cell reaction* 


Calculated open-circuit 
voltage at 25°C (V) 


Theoretical specific energy” (Whkg~!) 


Including O2 Excluding O2 
Li/O2 4Li +02 = 2Li20 2.91 5200 11,140 
Al/02 4Al +302 = 2Al203 2.13 4300 8,130 
Ca/O2 2Ca + O2 = 2CaO 3.12 2990 4,180 
Zn/O2 2Zn + O2 = 2ZnO 1.65 1090 1,350 


à The reduction of O to O?— usually occurs only in the presence of a catalyst; often the product is the peroxide, O77~. 
> Includes only the active materials. Since O2 does not have to be carried in the battery, values are given for the cases of including and excluding O2. The battery 


weight will increase once the discharge begins. 


were limited by oxygen diffusion and that dissolved oxygen 
diffusivity and concentration are of critical importance for the 
performance of a lithium/air or oxygen flooded-electrolyte cell. 
Oxygen concentration in an organic electrolyte is dependent on 
the electrolyte formulation and gaseous oxygen partial pressure. 
Ways to enhance the dissolved oxygen diffusivity and concen- 
tration in the organic electrolyte were by the use of different 
electrolyte formulations and by the increase of the gaseous oxy- 
gen partial pressure. 

Metal catalysts in the carbon/air cathodes were employed [3] 
to enhance the practical energy density of the lithium/air battery. 
The catalysts investigated were: manganese (Mn), cobalt (Co), 
ruthenium (Ru), platinum (Pt), silver (Ag), and a mixture of 
Co and Mn. The Mn catalyzed carbon cathodes showed the best 
performance followed by the Co catalyzed carbon cathodes. The 
specific discharge capacity of the Mn catalyzed carbon cathodes 
was determined to be 3137 mAh g7! of carbon at a discharge 


rate of 0.1 mA/CM com at a cut-off voltage of 1.5 V, whereas 


the Co catalyzed carbon cathode capacity was 2414 mAh g7! of 
carbon. 

Using in situ mass spectrometry [6], it has been shown that 
the electrochemical reaction given by Eq. (1) in a nonaqueous 
electrolyte, such as 1 M LiPF¢ in propylene carbonate (PC), is 
reversible and charge/discharge cycling can be carried out. The 
fundamental and practical issues are: the role of electrochemical 
manganese dioxide and other transition metal compounds in pro- 
moting the cathode electrode reaction, optimization of electrode 
porosity, structure, and composition of the cathode electrode. 
Prevention of water or carbon dioxide from entering the cell, 
when operated in air, is also a very important issue. 

A preliminary analysis of the lithium/air battery system was 
performed [7]. In this analysis, the reversible, individual elec- 
trode and overall cell reaction voltages for the formation of only 
Li2Oxs) and both Li2Oçs) and Li2O2çs) as the solid discharge 
product species in the air cathode reaction layer, amounts of the 
chemical species Li2O2(s), LizO(s), and Lis) as a function of cell 
discharge time, and the theoretical discharge capacity and ther- 
mal efficiency of the cell were determined. The total composite 
cathode volume for a desired maximum cell discharge, pore vol- 
ume, and volumes of the carbon and lithium peroxide materials 
were computed using developed analytical expressions. A math- 
ematical expression, based on the first law of thermodynamics, 
to compute heat exchange between the lithium cell or battery 
and its surroundings was also developed. 


Electrical load 


Carbon 
Cathode 


Separator 


Anode 
Protective 
Layer 


Fig. 1. Simple sketch of lithium/air system. 


The formulation presented in this paper is for a lithium/air 
cell, consisting of a thin lithium sheet anode, lithium phospho- 
rous oxynitride (LiPON) or glass anode protective layer, and a 
porous carbon oxygen/air cathode; see Fig. 1. The porous car- 
bon cathode is assumed to be made of carbon black, catalyzed 
by a metal such as manganese, polytetrafluoroethylene (PTFE) 
as the binder, and the pores are flooded with a lithium based 
electrolyte, for example, 1 M LiPF¢ in PC/1,2-dimethoxyethane 
(DME) in a mass ratio of 1:1 (1:1). 


2. Mathematical formulation 


Dependence of the lithium/air system performance on oxygen 
solubility, concentration, and its transport in organic electrolytes 
was experimentally shown [4,5]. A diffusion-limited model for 
oxygen transport and its conversion in the porous cathode reac- 
tion layer is presented below. 

Due to the low current densities, voltage losses associated 
with current collection, the lithium anode half-cell reaction, 
and lithium ion transport in the electrolyte are assumed to be 
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minimal. The flooded-electrolyte cell is assumed to be diffusion- 
limited with respect to dissolved oxygen. The porous carbon 
cathode is assumed to be catalyzed by a catalytic agent, for 
example, manganese (Mn). Pores in the carbon reaction layer 
are assumed to be full of the liquid phase electrolyte such as 1 M 
solution of LiPF¢ in an organic solvent mixture, e.g., PC/DME 
(1:1). The electrolyte has an initial concentration of dissolved 
oxygen. This initial oxygen concentration may be equal to zero 
(i.e. no oxygen present in the electrolyte) or equal to the solubil- 
ity limit (i.e. saturated concentration of oxygen at the prevailing 
temperature and pressure conditions). Furthermore, the con- 
dition of isothermality is assumed to prevail. During the cell 
discharge, oxygen diffuses through the liquid phase electrolyte 
in pores from the pure oxygen or air side of the cathode towards 
its side in contact with the lithium anode protective layer, for 
example, the LiPON or glass layer; assuming that the thickness 
of electrolyte separator layer between the protected anode and 
cathode is negligibly small. Oxygen is consumed by its reaction 
with the lithium ions in the presence of electrons at the active 
sites located at the carbon-electrolyte interface. It is assumed 
that the lithium ion concentration in the electrolyte solution is 
sufficiently high so that the electrochemical oxygen reduction 
can be safely assumed to be pseudo first-order with respect to 
the dissolved oxygen concentration. The developed oxygen con- 
tinuity equation, describing oxygen transport and its reduction 
in the porous carbon cathode, is given below 


0 (A | Ra (3) 


dz T OZ] 


: (eca) 
—(& = 
ot va 


where ¢ is the electrolyte filled pore volume fraction in the cath- 
ode reaction layer, ca = co, is the oxygen molar concentration 
in the liquid phase electrolyte, t is the cathode pore tortuosity 
factor, and D4 is the oxygen mass diffusivity in the liquid phase 
electrolyte. The oxygen conversion rate, Ra, at the active sites 
located at the pores surfaces per unit cathode volume is given 
by 


2e E 
Ra = | > | ksca = 2ks | —ca (4) 
rp rp 


where ks is the overall electrochemical oxygen reduction rate 
coefficient at a given cell operational temperature, and Fp is the 
mean pore radius of pores filled with the electrolyte solution. In 
the derivation of Eq. (3), it has been assumed that dca/dr = 0 
(i.e. pores are sufficiently narrow), and dca /30 = 0 (symmetry 
with respect to the angular coordinate 0 around the axis of a 
pore). Note that here 


[n]a F 
ks = kse exp RT 


_ —E, In |a° F 
= |ks,e,0 exp RT exp RT (5) 


where kg.¢9 is the intrinsic reaction rate coefficient associated 
with the nature of the surface reaction involving the reactant 
species, oxygen, lithium ions and electrons. Ep is the reaction 
activation energy, i.e., energy barrier to the surface reaction 
occurrence in the absence of the cathodic overvoltage, (when 


Le ee 
E 


z =0, c,=¢y, 


Fig. 2. Schematic of an individual pore in the cathode reaction layer. 


|n| =0.0 V), and &° is the charge transfer coefficient across 
the cathode—electrolyte interface and has a value of 1.0 if it is 
assumed that 2 electron transfer step is the rate-determining step 
in the overall intrinsic electrochemical oxygen reduction mech- 
anism. a° is equal to 0.5 for an assumed one electron transfer 
rate-determining step. 

It is here assumed that there are np open-ended pores, each 
of mean pore radius, Fp, corresponding to a unit geometric area 
perpendicular to the zı coordinate (see Fig. 2). The initial pore 
volume in the cathode reaction layer just before the start of the 
oxygen reduction process per unit geometric area perpendicular 
to the zı coordinate is given by 


(FÀ gTle)np = le£0 (6a) 


The pore volume in the cathode reaction layer during the oxygen 
reduction process is given by 


(nF te)Mp = kee (6b) 


Division of Eq. (6b) by Eq. (6a) leads to 


Py 
e = o( 7) (7) 


where £o and 7p,9 are the cathode void fraction and mean pore 
radius, respectively, just before the start of the oxygen reduction 
process, and, £ and 7p are the cathode void fraction and mean 
pore radius, respectively, at any time during the oxygen reduction 
process. It is assumed that the oxygen mass diffusivity in the 
liquid phase electrolyte and the pore tortuosity are constant. The 
following dimensionless quantities are now defined as: 


r=- (8a) 
Fp,0 
Cc 
caa = (8b) 
CA,0 
Z 
E= (8c) 


where C4 o is the oxygen molar concentration in the liquid phase 
electrolyte at the saturation level of dissolved oxygen or at the 
level of oxygen solubility in the electrolyte in contact with the 
pure oxygen reactant or air supplied to the cell cathode, le is 
the thickness of the porous cathode, and r, cag, and & are the 
dimensionless pore radius, oxygen concentration, and distance 
along the zı coordinate increasing towards the protective anode, 
respectively. Substituting for Ra and € from Eqs. (4) and (7), 
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respectively, into Eq. (3), and expressing Fp, ca, and z in terms 
of the dimensionless quantities, the following continuity equa- 
tion for oxygen in the dimensionless form is obtained 


(r ?caa) 9 (4 9cAa 
a) ; 9 
ay! 2E aE B(rca.a) (9) 
where 
, tDa 
A 9 
i (9a) 
and 
2k, t2 
ee (9b) 
Dafp,o 


The rate of decrease of pore volume per unit total porous cathode 
volume due to the formation of solid Li2zO2s) via the overall 
reaction, 2Li* + 2e~ + O2 = Li202(s), is given by 
ðe -RA Mti,0, 
ðt PLi2029 


(10) 


Substitution for Ra and € from Eqs. (4) and (7), respectively, 
into Eq. (10), and further simplification leads to 
OFp ksca Mj j,0, 


= (11) 


ot PLi2O25) 


Expressing Fp and ż in the dimensionless form, Eq. (11) becomes 


or 


af T YA (12) 
where 
"E (ksca, o MLiz03/PLi2029)) (12a) 


(Da?p,o/t2) 


To solve the coupled differential equations given by Eqs. (9) 
and (12) for ca q and r analytically or numerically, the following 
conditions in the dimensional form are used. 


e Initial conditions, t < 0; 0 < z1 < k: 


7p = Fp,0 (13a) 

ca = ch (13b) 
e Boundary conditions, t>0: 

zı =0, CA = CA0 (14a) 

ask, tL (14b) 


The above conditions are turned into the dimensionless form 
using the quantities given in Eqs. (8a)—(8c). 


e Initial conditions, ” <0;0<& <1: 


cee 


= (15a) 
rp,0 


r= 


ch 
Ad = n= (15b) 
CA,0 


e Boundary conditions, f > 0: 


o até=0 
wis Aa (16a) 
CA,0 
o até=1 
OCA,d 
— = 0 16b 
JE (16b) 


Oxygen molar flux into the electrolyte at any time during the 
cell discharge period, per unit geometric area perpendicular to 
z1, is given by 


i Daé OCA 
TE 


n 17 
T Oz] ny) 


21=0 


On substitution for £ from Eq. (7) into Eq. (17a) and the use of 
dimensionless quantities, one obtains 


Noo €o DACA,0 2 ðCA,d 
a tle oE é=0 


The current generation rate per unit geometric area perpendic- 
ular to zı in a differential thickness of dz,, via the occurrence 
of the electrochemical reaction: 2Li* +2e7 + O2 = Li20 xs), is 
given by 


(17b) 


digeom = 2FRa dz1 (18) 
Using the expression for Ra in Eq. (4) 
4Fk,e 
ieas | (18a) 
rp 


Substituting for (¢/7p) from Eq. (7) into Eq. (18a) and using the 
dimensionless quantities given in Eqs. (8a)—(8c), one obtains 


: Eole 
digeom = 4F(ksca,0) | =— } rea dé (19) 
rp,0 
Total geometric current density, corresponding to the oxygen 


reduction reaction occurring in the cathode reaction layer is then 
given by 


i geom 
lgeom = di geom 
0 


Eole =l 
= jræeao ( JI (rca,a) dé (20) 


Tp,0 


It should be noted that Eq. (20) gives oxygen reduction geometric 
current density at any time during the cell discharge. At the end 
of each time interval, the numerical solution of Eqs. (9) and (12) 
would provide data on the dimensionless oxygen concentration, 
cad, and the dimensionless average pore radius, r, as a function 
of the dimensionless distance, £. Using this data in Eq. (20), one 
can carry out the integration of Eq. (20) either numerically or 
graphically to determine the average geometric current density at 
the end of each time interval. If igeom is to be kept at a constant 
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Table 2 

Parameter values used in simulation 

Parameter Value(s) 

k 0.001, 0.010, and 0.070 cm 

CAO 6.44692E—07 g mol cm~? at 0.2 atm O2, 1.934076E—06 gmol cm~? at 0.6 atm O2, 3.22346E—06 g mol cm7? at 1 atm O% (all for 
1M LiPF¢ in PC/DME) 

£0 0.73 [4] 

Da 7.00E—6 cm? s~! [5] 

Fp,o 2.00E—7 cm (estimate) 

MW1Li02 45.8768 g (g mol)! 

PLi202 2.3 gem™ [7] 

T 1.17 (calculated from t = e5°>) 

igeom 5, 1, 0.5, 0.1, 0.05 mA cm~? 

T 298.15K 

a 0.0722 [5] for 1 M LiPF6 in PC/DME 

at 1 


level with respect to time, Eq. (20) suggests that the product 
|k, i. (rca,a) dé| must be kept constant. Because f, a (rca.a) dé 
is likely to decrease with time; therefore, ks must be increased to 
a level so that [ks i (rca,d) ag] = constant with respect to time. 


From Eq. (20), in conjunction with Eq. (5), igeom is expressed as 


: coke Inel ao F 
= |4F ks e exp | ———— 
Tgeom ( Po eao) ( s,e Xp ( RT 


(21a) 


The geometric current density at a time, t= t; or t’ = t4, is given 
by 


f _ (4Eeokcao) (, ox [nlna F 
geom, t = ZT s,e EXP RT 


é=1 
x | (rca,a) dé (21b) 
é=0 i 
1 
Division of Eq. (21a) by (21b) leads to 
igeom c c a F i. (rca,a) dE 
F = exp | (In| = [n ln) i (22) 
İgeom,tı { fo(rea.a) dé} 
ti 


In order that the geometric current density can be kept constant 
with respect to time, i.e., (igeom/igeom,t;) = 1, it requires that 


fo eae) dé on 


~ {fo(rea.ade} 


Eq. (23) suggests that if the average geometric current density 
at any time, t, is to be equal to that at a time, t=t,, where the 
oxygen reduction overvoltage is |n°|,,; the overvoltage |n°| 
must be adjusted and equal to that given by Eq. (23). Because 
the numerator of the right-hand side of Eq. (23) will decrease 
as the cell discharge time increases; therefore, the cathodic 
reaction overvoltage, |7°|, should increase with time to maintain 
the geometric current density at a desired constant level. 


(In| — | 


exp | RT 


ti 


3. Results and discussion 


Under numerically constant geometric current density and 
oxygen solubility in the organic electrolyte solution, Eqs. (9) 
and (10) were solved numerically using an explicit finite differ- 
ence technique. The relationship in Eq. (23) was used to keep 
igeom Constant by arbitrarily setting In°|r, = 0 and for the ini- 
tial conditions used knowing that the integral on the bottom of 
the right-hand side of Eq. (23) is equal to 1.0. Table 2 lists the 
parameters used in the simulation. 

Fig. 3 shows the dimensionless pore radius as a function of 
percent of total pore length and dimensionless time for a geo- 
metric current density of 0.5 mA cm~?, a cathode thickness of 
0.070 cm and an oxygen partial pressure of 1.0 atm. The highest 
dimensionless time represents the total discharge time for the 
cell. The simulated discharge would end when the relationship 
provided by Eq. (23) “blew-up”, that is, there was not enough 
dissolved oxygen to sustain the needed current. Fig. 4 has the 
same parameters as Fig. 3 except the geometric current den- 
sity is 0.10 mA cm™?. Note that at the lower current density, 
as shown by Fig. 4 in comparison to Fig. 3, a higher utiliza- 
tion of the electrode is realized as shown by the lower steepness 


Dimensionless Pore Radius 


0 10 20 30 40 50 60 70 80 90 100 
Distance Along Pore (%) 


Fig. 3. Dimensionless pore radius versus pore distance for igeom = 0.5 mA cm~?, 


le =0.07 cm and O7 partial pressure = 1.0 atm. 
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Dimensionless Pore Radius 
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Fig. 4. Dimensionless pore radius vs. pore distance for igeom =0.1 mA em7?, 
le =0.07 cm and Op partial pressure = 1.0 atm. 


of the dimensionless pore radius versus the dimensionless pore 
axial distance profile. Also note that the pore radius gradient 
in Fig. 4 is less than that for Fig. 3 due to the lower current 
density. 

With respect to Figs. 3 and 4, Figs. 5 and 6 show the cor- 
responding dimensionless oxygen concentration profiles in the 
pores. Note that the concentration gradient in Fig. 5 is higher 
than Fig. 6 due to the higher current density. Also, the oxygen 
concentration limiting the discharge time is clearly evident in 
Figs. 5 and 6. 

To determine the specific capacities (mAhg`~ of carbon) 
of the modeled cell, as a first approximation it was assumed 
that the solid volume fraction of the cathode, 0.27, was initially 
composed of only the carbon material with its density in the 
cathode being 2.2 g cm~?. Fig. 7 shows the cell specific capacity 
as a function of oxygen partial pressure and geometric current 
densities. This curve is very similar to measured experimental 
values (see Fig. 3 in [5]). 

Fig. 8 shows the specific capacity as a function of geometric 
current density and cathode thickness for this pore flooded- 
electrolyte design. Clearly, as the cathode thickness is reduced, 


1 


Dimensionlesss Oxygen 
Concentration 
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Fig. 5. Dimensionless oxygen concentration vs. pore distance for 
igeom = 0.5 mA cm~?, l = 0.07 cm and O3 partial pressure = 1.0 atm. 
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Fig. 6. Dimensionless oxygen concentration vs. pore distance for 
igeom = 0.1 mA cm~?, le =0.07 cm and O2 partial pressure = 1.0 atm. 


Specific Capacity (mAh/g) 


0.0 0.2 0.4 0.6 0.8 1.0 1.2 


Oxygen partial pressure (atm) 
Fig. 7. Lines are the calculated specific capacities (mAh g7! of carbon) for 
le = 0.07 cm as a function of oxygen partial pressure and current densities of (a) 
0.5, (b) 0.1, and (c) 0.05 mA cm7?. Symbols are the experimental values from 
Fig. 3 in [5]. 
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Fig. 8. Specific capacity (mAh g~! of carbon) as a function of current density 


for an oxygen partial pressure of 1.0 atm and /,’s of (a) 0.001, (b) 0.01, and (c) 
0.07 cm. 


S.S. Sandhu et al. / Journal of Power Sources 164 (2007) 365-371 371 


the cell performance as measured in terms of specific capacity 
for a given current density is greatly enhanced. 

It has clearly been shown above by simulation that Read’s 
[4,5] lithium/oxygen cells were flooded with electrolyte. The 
simulations represented by Fig. 7 are closely correlated with 
experimental results under the same conditions. Fig. 8 shows 
that thick cathode cell designs are not optimal and that a trade- 
off between starved and flooded electrolyte cell designs also 
needs to be performed. 


4. Conclusions 


A diffusion-limited model for lithium/air battery with an 
organic electrolyte was developed and the results from the model 
compare favorably with experimental results. The model shows, 
for a given cathode pore radius, the effect of oxygen partial 
pressure, current density, and cathode thickness on the specific 
capacity of the battery. 
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